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Abstract 

In this work, we find the expressions of continuum HQET four-fermion operators in 
terms of lattice operators in perturbation theory. To do so, we calculate the one-loop 
continuum-lattice HQET matching for the complete basis of IS.B = 2 and /S.B = 
operators (excluding penguin diagrams), extending and completing previous studies. 
We have also corrected some errors in previous evaluations of the matching for the 
operator Oll- Our results are relevant to the lattice computation of the values of 
unknown hadronic matrix elements which enter in many very important theoretical 
predictions in i?-meson phenomenology: B^-B^ mixing, tb and tb^ lifetimes, SUSY 
effects in AS = 2 transitions and the Bg width difference AF^^. We have reanalyzed 
our lattice data for the Bb parameter of the B^-B^ mixing on 600 lattices of size 
24^ X 40 at /? = 6.0 computed with the SW-Clover and HQET lattice actions. We 
have used the correct lattice-continuum matching factors and boosted perturbation 
theory with tadpole improved heavy-light operators to reduce the systematic error in 
the evaluation of the renormalization constants. Our best estimate of the renormal- 
ization scale independent S-parameter is Bb = 1-29 it 0.08 it 0.06, where the first 
error is statistical and the second is systematic coming from the uncertainty in the 
determination of the renormalization constants. Our result is in good agreement with 
previous results obtained by extrapolating Wilson data. As a byproduct, we also ob- 
tain the complete one-loop anomalous dimension matrix for four-fermion operators in 
the HQET. 
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1 Introduction and Motivation. 



B-hadron decays are a very important source of information on the physics of the 
standard model (SM) and beyond. In many important cases, long distance strong con- 
tributions to these processes can be separated into matrix elements of local operators. 
Lattice QCD can then be used to compute these non-perturbative parameters from first 
principles. A list of some four-fermion operators relevant to i?-meson phenomenology 
is the following: 

AB = 2 operators 

Oll{rr) = 57^(l=F75)g&7M(l=F75)g 

OIl{rr) = ^(lT75)gft(lT75)g 

Olr{rl) = &7^(l=F75)g&7M(l±75)g 

OIr^rl) = 6(lT75)g6(l±75)g (1) 
AB = operators 

Qll(rr) = bj^ {l=f j5)q q'jf, {I =fj5)b 

Qll(rr) = ^(lT75)gg(lT75)& 

Qlr{rl) = b7^(l=F 75) 9^7/^(1 ±75)^ 

QlRiRL) = 6(lT75)gg(l±75)& (2) 

where q denotes a light quark n, d or s. The corresponding operators with a (the 
generator of the SU(3) group) insertion, denoted by OtiL ■ ■ - etc, will also be consid- 
ered. 

As is well known, Oll determines the theoretical prediction of the B^-B^ mixing 
in the SM. Moreover, Oll together with Of^p^^j^^, Orr and Of^, parameterize SUSY 
effects in AB = 2 transitions Oll and Of^^^ determine the Bg width difference 
AFb, [||. The AB = operators Qll and Qlr, parameterize spectator effects in the 
tb and tb^ lifetimes For the four-fermion operators with a t"" insertion: Ot^j^{RR) 
and Oi|^ contribute to SUSY effects in AB = 2 transitions and QIll and Qt^R 
to spectator effects in the B-meson lifetimes |^]. 

Our aim is to non-perturbatively evaluate the matrix elements between S-meson 
states of all the operators in eqs. (I) and (D, using lattice simulations of the 6-quark in 
the HQET. These non perturbative parameters can also be measured in principle with 
propagating h quarks by simulating in the charm mass region and extrapolating to the 
h quark mass. The procedure to perform the transition from QCD to lattice HQET, 
a combination of analytic and numerical calculations, can be split into the following 
four steps: 

First Step: the continuum QCD - HQET matching. 

The QCD operators are expressed as linear combinations of HQET ones in the 
continuum at a given high scale, say, ^ = rrii,. 
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Second Step: running down to fi = a in the HQET. 

The HQET operators obtained in step 1 at the high scale ji = rub are evolved 
down to a lower scale /x = a~^, appropriate for lattice simulations, using the 
HQET NLO renormalization group equations. 

Third Step: continuum-lattice HQET matching. 

Having obtained the continuum HQET operators at the scale /i = , they are 
expressed as a linear combination of lattice HQET operators at this scale. The 
procedure is very similar to step 1: two amplitudes are matched at one-loop order, 
one in the continuum HQET and the other in the lattice HQET. 

Fourth Step: lattice computation of the matrix elements. 

The matrix elements of the lattice HQET operators in Step 3, are measured 
by Monte Carlo numerical simulations on the lattice. Using these values and the 
chain of matching equations in steps 1 to 3, the matrix elements of the continuum 
QCD operators relevant to phenomenology can be determined. 

In this work, we deal with step 3, i.e. we find the expressions of the continuum 
HQET operators in terms of lattice HQET ones using continuum and lattice perturba- 
tion theory at one loop. The calculation of steps 1 and 2, and the numerical simulation 
in step 4 are in progress and will be published elsewhere. 

We would like to stress that a non perturbative determination of the lattice renor- 
malization constants would be preferable because at the values of the coupling at which 
the simulations are performed, the one-loop perturbative corrections are not small, due 
to the appearance of tadpole diagrams. In principle, the Martinelli's et al |0] method 
for non perturbative renormalization may be utilized to avoid the problems with per- 
turbation theory. Unfortunately, however, the application of the former to the HQET 
is not straightforward [|[. The key observation here is that HQET lattice simulations 
suggest that the finite heavy quark propagator in the Landau gauge has the form |^ 

S{x,t) = 6{x)e{t)A{t)e-^^-^'^^' (3) 

where A(t) is a smooth unknown function of t, A is a constant and 6fh is the residual 
mass (also called the HQET mass counter-term). The latter is needed to remove the 
linear divergence in A. On the one hand, the difference ms = A — (5m should be small 
in order to reduce 0{msa) terms. On the other hand, however, a too small ms implies 
that the heavy quark propagator goes to a constant. This fact makes its numerical 
Fourier transform, a necessary ingredient in the Martinelli's et al method, very prob- 
lematic due to the appearance of huge finite time effects of 0(6"™^"^), where T is 
the lattice time length. Therefore, the non perturbative calculation of renormalization 
constants in the HQET is difficult and requires more careful studies. 

The paper is organized as follows. We begin by introducing the lagrangian of 
the HQET in the continuum and its discretized version which we will use in lattice 
simulations. The corresponding Feynman rules are derived for both heavy quarks 
and antiquarks. In section 3, we give the lattice lagrangian for light quarks and the 
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corresponding Feynman rules. The strategy of the calculation is presented in section 
4. In this section, our renormalization prescription is discussed in detail. Moreover, 
we study the subtle point of the presence of evanescent operators and their effect on 
the matching and higher-order calculations. Then, we present the calculation details 
of the Feynman diagrams in the continuum and on the lattice HQET using several 
renormalization schemes. In section 6, our results for the one-loop continuum-lattice 
matching are given and compared with previous partial calculations. We also discuss 
the effects of a change of the continuum renormalization scheme. Then we reanalyze the 
mixing using the correct matching factors. Finally, we present our conclusions 
in section 8. In addition, we include one appendix where the relevant formula used in 
our calculation are collected. 



2 Feynman rules for the continuum and discre- 
tized HQET lagrangians. 

The Minkowski continuum lagrangian for a field h(x) which annihilates a static heavy 
quark and for a field h{x) which annihilates a static heavy antiquark is 

^hqet{x) = h\x)iDoh{x) + h{x)iDoh^{x) (4) 

where zD^ is the covariant derivative idfj, + gt"" A'^{x) with the gluon field and 

= the generators of the SU(3) color group normalized by tr{t"-t^) = 5"'^/2. A 
very important property of the HQET fields h and h is the fact that joh = h and 
/i 7o = — h. We will make extensive use of these properties in our calculation. 

On the lattice, in Euclidean space, we have first to choose a discretization for the 
continuum covariant derivative D^, which has the correct continuum limit. Obviously 
there are many equivalent possibilities. The point is that the discretization we use in 
the perturbative calculation of the lattice operators must be the same as the one with 
which we perform the numerical simulation. Otherwise, the mismatch will generate 
spurious uncontrolled finite contributions to the final results. Our convention is to 
discretize backward in time, as first proposed by Eichten 0], so that the lattice HQET 
action is 

gHQET = a^Y,{ hHn) [h{n) - C/](n - 4) h{n - 4) 



+ 



h{n) [c/4(n)/i^(n + 4) - /i^(n)]} (5) 



where n denotes a lattice site, 4 is the unit vector in the Euclidean time direction and 
U is defined in eq.@ bellow. 

From eqs.(Q) and (||), it is easy to find the Feynman rules which are collected in 
Table |. Some remarks are in order here. To perform a weak-coupling expansion in 
terms of g of the lattice lagrangian, we have parameterized the link variables as 

;7t(x) = 6-*"^^^^*" (6) 
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Quark 


Minkowski HQET 


Lattice HQET 


Propagator 

quark-gluon vertex 
seagull vertex 


i 

+ ie 
does not exist 


a 

1 + e - e-*P4a 


Antiquark 


Minkowski HQET 


Lattice HQET 


Propagator 

quark-gluon vertex 
seagull vertex 


i 

+ ie 
does not exist 


a 

1 + e - e-*P'*'* 
+ igt'^,p^^.Ae-'P> 



Table 1: Feynman rules for the continuum Minkowski HQET and the discretized Euclidean 
HQET. Our conventions are discussed in the text. 

where A^j^^x) is the gluon field and a is the lattice spacing. In Table [l|, the propaga- 
tor corresponds to a quark (antiquark) of momentum and a color conservation is 
understood. The quark-gluon vertex corresponds to an incoming quark (antiquark) of 
color a and momentum outgoing quark (antiquark) of color (3 and momentum p' 
and incoming gluon of color index a and Lorentz index /x. The lattice seagull vertex, 
which will give rise to tadpole diagrams, is defined as the quark-gluon vertex except 
for an additional incoming gluon with color index h and Lorentz index u. 



3 The light— quark Wilson and Clover actions. 

The operators in eqs.(|l|) and (Q) contain both a heavy and a light quark. The dynamics 
of the former is governed by the HQET lagrangian described in section 2. In this 
section, we consider the dynamics of the light degrees of freedom. 

In the continuum, we use the standard QCD lagrangian for light quarks, which 
Feynman rules can be found in many text books (see, for example, ref.|^]). 

On the lattice, and in our numerical simulations, light quark propagation is de- 
scribed through the Wilson action 

nil ^ ^ 

Ar 1 

(7) 



+ V5(n + /i)(r + 7^)[/t(n)V(n)] } + 5] 



4r 

V'(n) (mo H ) i>{n) 

a 



where is the standard Wilson plaquette action r is the Wilson parameter 

and mo is the bare light fermion mass. 
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An important source of systematic error in lattice simulations is the finiteness of the 
lattice spacing, a. It is well-known that hadronic matrix elements computed with the 
Wilson action differ from the corresponding continuum ones by terms of 0{a). Some 
methods have been devised to reduce the cut-off dependence of the lattice hadronic 
matrix elements involving Wilson fermions. In this paper, we will apply the method of 
improvement first proposed by B. Sheikholeslami and R. Wohlert (SW) long ago jlO[ |, 
which has been deeply studied in ref.[|l^]. The main idea of this technique is to use the 
freedom we have in the definition of the discretized quark action to eliminate all terms 
of 0(a) from the lattice matrix elements, by adding some irrelevant 0{a) terms to the 
action (^) and redefining the lattice fields ip. 

The first recipe is to perform all calculations with the so-called SW-Clover action, 
given by 

. . r nr — 

(8) 



QjT — 



n,ij,,u 

where P^u{n) is the discretized field strength tensor defined by 

1=1 ^ 

and the sum is over the four plaquettes in the ^-v plane, stemming from the point n 
and taken counterclockwise. 

The second recipe is the rule for constructing improved operators: the quark fields 
in lattice operators must be rotated 



(1 - z) mo 



(10) 



where ^ is the symmetric lattice covariant derivative and z is an arbitrary real number. 

It can be demonstrated that on-shell lattice matrix elements of improved operators 
computed with the SW-Clover action do not contain 0{{g'^)"'a log" a) terms for any 
value of z. We will take z = 1 so that the rotated operator do not depend on rriQ. The 
reason for this choice is that in this case perturbative calculations are more accurate 



0. 



The same method can be applied to the Eichten's action for heavy quarks in eq.(|5|). 
In ref.|T^, it is shown, however, that for 0{a) improvement of on-shell matrix ele- 
ments, no modification of the static quark propagator is needed. Therefore, in order to 
eliminate 0{a) terms from on-shell matrix elements of heavy-light operators, we have 
to perform the calculation with the SW-Clover action and rotate the light quark fields 
only. For heavy quarks, the usual unimproved Eichten's action dH) can be used. 

For practical reasons, instead of rotating the fields in the numerical simulation, a ro- 



tation of the light-quark propagator is done [13]. This technique is equivalent to a new 



definition of the lattice quark fields and automatically improves the operators. Since 
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sm(p^a)) + (mo + 2 £ sin2((p^a)/2)) ] 


Wilson quark-gluon vert. 




7p cos((p + p')pf) - sin((p + y)p§) 


e ^p'p 2 e*^'' 2 


Impr. quark-gluon vert. 




E,. ( o-pu sm{{p - p')„a) ) cos((p - p')pf ) 


e~*^p'2 e*Pp2 


Impr. quark-quark vert. 




quark-gluon-quark vert. 





Table 2: Feynman rules for the lattice Wilson and SW-Clover actions for light quarks. 



this method is much easier to be implemented on the computer than the improvement 
with rotated fields, we will apply it in this paper. 

To be definite, let us consider a generic AB = 2 effective operator with arbitrary 
Ti and r2 dirac matrices 



0(n) = (^h\n)Ti 



(11) 



The improved operator is obtained by rotating the light-quark fields according to 
eg. (p!o[) , keeping terms of order O(a^) for the second method of improvement |13], 



+ 



h\n)Tiij{n)) f /i(n) Ta ^ V(n) 



h^{n)Ti]^ij{n)j f /i(n) V(^) 



(12) 



Similarly, we can improve the AB = operators. 

From eqs.(^), (^) and (|l2|), we can obtain the Feynman rules which are given in 
Table ||. In this Table, the propagator corresponds to a quark of momentum p and bare 
mass mo, and color conservation is understood. Both the Wilson and the improved 
quark-gluon vertex, the latter coming from the SW term in action (|8|), correspond 
to an incoming quark of color a and momentum p, outgoing quark of color /? and 
momentum p' and incoming gluon of color index a and Lorentz index p. The other 
two vertices in Table 2 come from the rotation of the four-quark operators in eq.(|T2|). 
In the improved quark-quark vertex we have an incoming or outgoing light-quark of 
color a and momentum p. In the quark-gluon-quark vertex, we have in addition an 
incoming gluon with momentum q, color index a and Lorentz index p and the sign of 
the momentum p in the exponentials is + for an incoming light quark and — for an 
outgoing one. For the gluon, we use the standard Wilson action 0]. 
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4 The strategy: the subtle point of the renor- 
malization scheme dependence. 

Our aim is to express the continuum HQET operators at the scale fi = as a hnear 
combination of lattice HQET operators at this scale in such a way that all continuum 
matrix elements be equal to their corresponding lattice counterparts up to one-loop 
order in perturbation theory. 

The procedure is well known. The relation between the continuum and lattice 
operators to one-loop can be written as 

or = Y.i^^3 + (13) 

By sandwiching eq.([l^) between appropriate, but arbitrary, initial and final states, one 
can obtain the renormalization constants Z. In fact, some amplitudes are evaluated 
up to one-loop in perturbation theory both in the continuum and the lattice theory. 
After renormalization one gets 

2 

j 

where (0) denotes the bare matrix element. Thus, the constants Z are obtained by 
simply subtracting at the scale given, Zij = cf-^ ~ c^if\ this is the matching procedure. 

Note that the amplitudes can be regularized and renormalized using, in general, 
different schemes in both theories. Some subtle points to be stressed are the following: 

1. the matching must be calculated with the same action as for the numerical sim- 
ulation because the matching coefficients depend on the action to be simulated. 

2. due to the breaking of chiral symmetry induced by the Wilson term for light 
quarks, the original operators can mix with lattice operators with different chi- 
rality. 

3. the matching constants Z's depend on the renormalization procedure chosen to 
define the operators in the continuum. This dependence will cancel with the 
remaining scheme dependence of the running in step 2 of the Introduction. This 
subtle point will be studied in detail in a forthcoming paper. 

Therefore, it is very important to clearly define the regularization and renormalization 
schemes we utilize in our calculation. Many prescriptions can be chosen to regular- 
ize the corresponding Feynman integrals: NDR (Naive Dimensional Regularization), 
DRED (Dimensional Reduction), HV (the 't Hooft-Veltman scheme) ...etc. Differ- 
ent regularization schemes will give rise to different finite parts in the amplitudes and 
hence to different matching coefficients. Having selected the regularization procedure, 
we have to choose the renormalization scheme to eliminate ultraviolet divergences, for 
example MOM (momentum subtraction), MS . . . etc. 
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Apparently, the continuum operators are now completely defined. A more careful 
analysis demonstrates, however, that this is not the case. In fact, there is another 
prescription to be set, namely, how to deal with the evanescent operators which un- 
avoidably are generated We will explicitly show that different treatments of the 
evanescent operators will give rise to different finite contributions to the amplitude even 
when we return to four dimensions. The mechanism is clear: evanescent operators are 
operators which vanish in four dimensions, in other words, they are operators of order 
0(n — 4), with n = 4 + 2e. When an evanescent operator is produced in a divergent 
diagram, terms of the form 1/e x evanescent arise. Therefore, 0{e^) contributions are 
generated and their value depend on the prescription we use to subtract the evanescent 
operators. 

To sum, it is not sufficient to specify the scheme we use to renormalize our Feynman 
diagrams, we have also to define the prescription to eliminate the evanescent operators 
Ip!^ . This fact results in the existence of families of MS schemes which members differ 
in the treatment applied to the evanescent operators: we can, for instance, simply 
subtract them together with the divergent pole or project the amplitudes onto some 
convenient, but arbitrary, operator basis using some well-defined projectors. Note 
that in the second case, the projectors can be chosen in many different ways, each of 
which will define a different member in the corresponding scheme family. Of course, 
we can use our favourite prescription. It is immaterial with which one we perform the 
calculation. At the end, when the continuum QCD operator will be constructed in 
terms of HQET lattice operators, all dependence on the intermediate renormalization 
scheme used to define continuum HQET operators will drop and the final result will 
depend on the definition of the QCD operators only. The key point here is to be 
consistent during the whole procedure so that no spurious scheme dependence creep 
into our calculation. 

5 The calculation. 

We compute the matrix elements of the AB = 2 operators in eq.(|l|) between an initial 
B^{bq) state and a B^{qb) final state. For the AB = operators in eq.(||) we calculate 
the matrix elements between an initial and a final B^{qb) state neglecting penguin dia- 
grams 1^]. We set the light-quark mass to zero and take our momentum configuration 
to be vanishing incoming and outgoing momenta of both light and heavy quarks. To 
eliminate the infrared divergences which appear at zero external momenta, we give the 
gluon a mass, A. This procedure is justified at one loop order. 

On the lattice, the Feynman diagrams to be evaluated are depicted in fig. 1. The 
box represents an insertion of the four-fermion operator. The cross vertex stands for 
the contribution coming from the rotated part of the improved operator (see the entries 
for the improved quark-quark and quark-gluon-quark vertices in Table Q). Note that 
each light quark-gluon vertices in all diagrams consists of the Wilson vertex plus the 
improved contribution coming from the Clover action (see entries for the Wilson and 
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Improved quark-gluon vertices in Table ^). Ultraviolet divergences are regulated by 
the lattice cut-ofF, a. 

In practice, we project the lattice amplitudes onto the Dirac basis B = {7zp, 7^±, 
7^, 7^±, o"=p }, where 7::^ denotes the Dirac structure 7'^(1 =F 75) ® 7^(1 =F 75)) 7q:± is 
7''(1=F75)«'7/.(1±75), 7^ is (1 =F 75) (1 =F 75), 7^± is (1 =F 75) ® (1 ± 75) and finally 
azf is =F 75) ® i^/ii/Cl ^75)- Notice also that a^:± vanishes in four dimensions. The 

projectors are defined by (the sign is + in Euclidean space and — in Minkowski space) 



p^^(ri®r2) 
p^^(ri®r2) 



where n 



Tr 
Tr 
Tr 
Tr 
Tr 
Tr 
Tr 



ri7'^(l ± 75)r27M(l ± 75)]/(32n(2-n)) 
ri(l T 75)r2(l ± 75)]/(32n) 

ri ( 1 T 75 ) r2 ( 1 T 75 ) ] i-n" + 9n - 16)/(256 (n - 2)) 
Tia'^'il T 75)r2fT^.(l T 75)]/(256(2 -n)) 
ri7'^(l ± 75)r27^(l T 75)]/(32n) 
Ti a'^" ( 1 T 75 ) r2 ( 1 T 75 ) ] /(256n (n - 1) (2 - n)) 
ri(l T 75)r2(l T 75)]/(256(2 -n)) (15) 



4 (n 



4 + 2e) for a four (n) dimensional gamma algebra and ri^2 
are arbitrary Dirac matrices. On the lattice, 7^ are the Euclidean gamma matri- 
ces satisfying the anticommutation relation {7^,7,/} = 2 5^1, in four dimensions and 
C/ii/ = 1/2 [7/i, 7i/]- With our definition, all operators in the basis B project back onto 
themselves. 

Since HQET fields h satisfy ^qUh = Uh and 'joVh = —Vh, where Uh and Vh are the 
spinors for a heavy quark and antiquark, we can reduce the number of independent 
lattice amplitudes through the following Euclidean-space relationships 

[uhcrf"' {I =F l5)vgVhaf,„{l =F 75) ^^<?] = 

-4 {[Uh{l =F l5)VqVh{l =F 75) Uq] + [Uhl^C^ T I5) Vq Vh IfL {'^ T 75) ^^g]) 
[Uh Cr'"' (1 =F 75) Vq Vq G (1 =F 75) Uh] = (16) 
+4 {[Uh{l =F l5)VqVq{l =F Uh] - [Uhl^ {I =F 75)^^g^^q7^ (1 ± ^5) Uh]) 

where Uq and Vq are the light quark and antiquark spinors. The first relation can be 
used to reduce Ai? = 2 amplitudes and the second for AS = ones. 

Moreover, since Fierz transformations are well defined in four dimensions, we have 
for A.B = 2 operators 



^[7''(1 T 75) ® 7m(1 T 75)] 
T[-f^{l T 75) ® 7m(1 ± 75)] 
■^[(1 T 75) €5 (1 T 75)] 



+ 



- [7^ (1 T 75) ® 7m (1 T 75)] 
2 [(1 ± 75) (1 T 75)] 
[(1 T 75) ® (1 T 75)] 

\ b'' (1 T 75) 7a. (1 T 75)] 



(17) 



where J- denotes the Fierz transformation and we have omitted for simplicity the 
external spinors. 
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In the continuum, only diagrams Di to Dq contribute to the amphtude, for neither 
improved nor tadpole vertices appear (see Table 0). Ultraviolet divergences will be reg- 
ulated using dimensional regularization. Unlike the lattice case, evanescent operators 
do arise in the continuum so that we have to give a prescription to deal with them. 

If the gamma algebra is in n dimensions, as in NDR, the relations in eq.(^) are not 
valid because Fierz transformations are not well defined. Moreover, since the "magic" 
formulae for reducing the product of three gamma matrices is correct only in four 
dimensions, eq.(p^) holds up to evanescent operator contributions. For instance, the 
continuation of the first relation in eq. (|lq) to n dimensions can be defined by 



(1 T 75) ® (1 T 75)] = 

/(n) ([7^ (1 T 75) ® 7/. (1 T 75)] + [(1 T 75) ® (1 T 75)]) + c^f (18) 

where f{n) depends on the choice of the evanescent operator uj. A naive prescription 
is to impose that f{n) = 4 as in four dimensions (notice that there is a change of sign 
between Minkowski and Euclidean spaces). However, we can equally well take /(n) = n 
by redefining cTf. A very interesting choice is /(n) = 4(1 + 1/6 e). This prescription 
is useful because leads to results which are Fierz symmetric. Notice that all three 
prescriptions are valid and differ only in the definition of the evanescent operator cTf. 

If we use DRED, where the positions and momenta are continued to n dimen- 
sions whereas all other tensors, in particular the gamma algebra, are leaving in four 
dimensions |15|, eqs.(|l^) and ( |T7|) hold but new independent operators arise |16| 



7^ = [7^ (1 T 75) 7" (1 T 75)] 5^"^ = 4 7t + 7^ 
7t± = [7^ (1 T 75) 7" (1 ± 75)] g^^J = J 7t± + lU 



= [a^^'^ (1 T 75) ® <y^^ (1 T 75)] gf, g^^ ^ J a± + 4 (19) 

where g"^ and gl^v are the metric tensors in four and n dimensions respectively. The 

Dirac structures 7^, 7zp± and are defined similarly to those with a bar except for 

the fact that Lorentz indices are contracted in four dimensions only. We have also 

f f f 

defined three evanescent operators 7^, 7^_|_ and which will prove to be useful in the 



calculation. Notice that the operator a±^, with the same notation as in eq.(19), is an 
evanescent operator because a±^: vanishes in four dimensions. 

In the scheme proposed by 't Hooft and Veltman (HV) |jl^, 75 anticommutes with 
the gamma matrices in four dimensions but commutes with those in n — 4 dimensions. 
In addition, the rule to consistently define the coupling to chiral fields in the Stan- 
dard Model implies that Dirac matrices in the four-fermion operators must be in four 
dimensions. 

In order to show the dependence of the matching on the prescription chosen to 
define the continuum operators, we will perform the calculation in several schemes, 
namely. 
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NDR^: NDR-MS, where we project onto the basis B in n dimensions with eq.(15). 
The HQET magic relation for AB = 2 operators (|l^), used to reduce the number 
of independent amphtudes, is taken with f{n) = 4(1 + ^e), where ^ is an 
arbitrary real number. A similar prescription is used for AB = operators. 



DRED I: DRED-MS', where we subtract with the basis B plus the evanescent oper- 
ators {7t - 7t' 7=f± - 7t±' - ^T' ^T±}- 



DRED II: DRED-MS", where we subtract with the basis B plus the evanescent op- 
erators {7^, 7^±, ai^, a^± }. 



HV: HV-MS, where we project onto the basis B in four dimensions with eq.(|l5|). 

Notice that DRED II is equivalent to eliminate the evanescent operators by simply 
projecting them out because -Pi (7/) = (7^ being any of the evanescent operators in 
the basis) for all projectors Pi in four dimensions in eq.(ll5|). 



6 The Results. 

The continuum HQET operators Oll, Olr, Of^ and Of^ can be written in terms of 
lattice operators as 
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where the analytical expressions for the lattice constants, D, can be found in the 
appendix A. Their numerical values for the Wilson parameter r = 1.0 are given in 
Table ^. In obtaining these numbers, we have used a reduced value of the heavy-quark 
self-energy constant e = 4.53 (to be compared with the non reduced one e = 24.48) 
which is consistent with fitting the time dependence of the Monte Carlo data through 
Ae^\ Notice that we have explicitly separated the contributions coming from the 
Wilson action (denoted by Dll ■ ■ ■ etc) , the Clover action with improved operators 
not including O(a^) terms (denoted by ...etc) and including O(a^) corrections 
(denoted by Z)£i . . . etc). The new lattice operators Oat, Om, Op and Oq are defined 
by 
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As can be seen from eqs.(|20|) to ([23|), the lattice transcriptions of the continuum 
operators depend on the renormalization scheme used to define the latter, including 
the prescription to deal with the evanescent operators. If fact, z to z'" parameterize 
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Wilson 
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Impr. II 
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Wilson 


Impr. I 


Impr. II 
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Table 3: Numerical values of the constants D for the Wilson parameter r = 1.0. 
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NDR^ 
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DRED I 














DRED II 





1 








HV 











1 



Table 4: Values of the scheme dependent parameters z's. 



the scheme dependence in the continuum and are defined in Table A interesting par- 
ticular value of ^ is ^ = 1/3 which corresponds to a NDR scheme which preserves Fierz 
symmetry, i.e. renormalization and Fierz transformations commutate in the continuum 
scheme NDR1/3 at one-loop order. 

Notice that the results in eqs.(^) to (23) are expressed in terms of two coupling 
constants, namely, the continuum, Qs, and the lattice, giat, ones. At one loop, how- 
ever, we can consistently identify these two coupling constants because the difference is 
0{a1). In this case, the logarithms combine giving log(a^ fj?) so that the infrared regu- 
lator A disappears, as it should be. Several recipes for improving the convergence of the 
lattice perturbative series by choosing a convenient perturbative expansion parameter 
have been proposed [19|. The choice of the optimal coupling constant is discussed for 
the particular case of the operator Oll in the next section. 

The lattice counterpart of the operator Oll is already known. It was determined 
by Flynn et al |21| for the Wilson action and by Borrelli and Pittori [12| for the SW- 
Clover action. We have extended the calculation to all the HQET operators in eq.([l|). 

but disagree with ref.|l^ 



Our results for Oll agree with ref.[21| but disagree with ref.|I^ in the sign of the 
contribution proportional to Oj^|j coming from the Feynman diagram Dq in fig. 1. 
This error propagates to the constant -D/j^ which correct value is given in Table | (to 
be compared with the value quoted in ref.[12|, Dj^j^ = —2.58.) 

The operators with a insertion, denoted by Oti, can be readily obtained from the 
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corresponding expressions in eqs.(^) to (^) by simply replacing the lattice operators 
Oi with the operators Otj and the operators On, Om, Op and Oq, defined in eq.(24), 
with On', Om', Opi and Oqi given by 
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For the continuum HQET Ai? = operators, the results we have obtained are the 
following 
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where the analytical expressions for the lattice constants E can be found in the ap- 
pendix A. Their numerical values for the Wilson parameter r = 1.0 and a reduced 
heavy-quark self-energy, are given in Table ^. Notice that, as before, we have ex- 
plicitly separated the contributions coming from the Wilson action (denoted by Ell 
...etc), the Clover action with improved operators not including O(a^) terms (de- 
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Table 5: Numerical values of the constants E for the Wilson parameter r = 1.0. 



noted by . . . etc) and including O(a^) corrections (denoted by Ej^^^ . . . etc). The 
new lattice operators Qr, Qs, Qt and Qu are defined by 
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The operators QtR^s,T,Uj with a insertion, are defined as those in eq.(|30D simply 
replacing Q with Qt. 

The expressions of the AB = continuum HQET operators with a can be 
obtained from eqs.(p6|) to ( |29| ) through the following procedure. Let Qt be a generic 
continuum Q-operator with a t° insertion. In terms of lattice operators, it can be 
written as 

Qt = Y, Eti Qtf + Y.Ei 

i i 

where Fti and Fi are related to the corresponding constants of Q as follows 



Ft 

Fi = 2/9 



Fti = -1/2 Eti 
Fi = l/QEi 



Hi =R, S, T, U (31) 
As a byproduct of our calculation, we have also obtained the anomalous dimension 
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matrix, 7, at one loop, defined by 
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(32) 



for both the AB = 2 and AB = operators. In the continuum dimensional regular- 
ization 7(0) is simply twice the residue of the 1/e pole and on the lattice it is twice the 
coefficient of the log(A^ a^). Therefore, we get for the AB = 2 operators. 



7 



(0) 



(1 + i) 







-37V + 4 + ^ 









3/2 






6 

-3(iV- 



2^- 

N' 



in the basis {Oll, Of^^, Olr, O^^} and where is the number of colors. 
For the AB = operators, we obtain the submatrices 
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(34) 



in the basis {Qll, QtiL-. Qlw Q^lr}^ and 
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(35) 



in the basis {Qlr, QtiR, Qll^ Qthh)- 



7 Reanalysis of the B -B mixing 

As we said in the previous section, the value of the lattice constant -Djj/j quoted in 
ref.[^], is incorrect. This constant enters the continuum-lattice matching of the QCD 
operator Oll, which matrix element between B^ and B^ states determines the the- 
oretical prediction of the B^-B^ mixing through the so-called B-parameter Bb ||2^ . 
In this section, we re-calculate the relevant renormalization constants and give the 
correct value of the B-parameter. 
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The B-parameter Bb is the ratio of the matrix element of the operator Oll to its 
vacuum insertion approximation 

, f I O^^O-) I (36) 
I I (0 I A„ I B»> I' 

The continuum QCD-lattice HQET matching for the operator Oll and the heavy-light 
axial current = h^^^^q can be written as 

Ull = ^Ol Oll + ^Or + ^Om On + ^Os Oll 

= ZAAf (37) 



where the Z^'s are matching (renormalization) factors which explicit expressions up to 
one-loop order are obtained by combining the results of steps 1 and 2 in the Intro- 
duction (we refer the reader to ref.|^2| and references therein, where all the relevant 
formula are collected 0) and eq.(^) which corresponds to step 3 in the Introduction. 
The values of the constants D's which enter eq.(|20|) are given in Table ^. Notice 
that in all previous Clover HQET computations Bb, the incorrect value of the con- 
stant Dji = Drr + D^f^^ + -D/yj = —5.4 was used, to be compared to the correct one 
Dli = —3.23. The physical value of the B-parameter is then given by 

i=LL,RR,N,S g I \U I ^4 V«J I / I i 

The numerical determination of the relevant lattice matrix elements, step 4 in the 



Introduction, is also discussed in detail in ref.|22|. We will take the values quoted in 
ref . p^] as input for our reanalysis. 

We now turn to the delicate point of the improvement of our lattice perturbation 
theory results. We start by discussing the choice of the expansion parameter. Many dif- 
ferent, but equivalent at one-loop level, definitions of the lattice coupling constant can 
be used in the perturbative calculation of the lattice-continuum matching. At scales 
~ 2 — 4 GeV of our numerical simulations, we expect small non-perturbative effects 
in the renormalization constants because of asymptotic freedom. However, tadpole di- 
agrams, which appear in lattice perturbation theory, give rise to large corrections and 
then to large uncertainties in the matching procedure. In refs.|T^, p^ ], some recipes, 
usually refer to as Boosted Perturbation Theory (BPT) , for choosing an optimal lattice 
coupling constant which would absorb many of the unwanted gluonic tadpole contribu- 
tions are suggested. The claim is that the perturbative estimate of the renormalization 
constants obtained from BPT at low orders, is closer to the non-perturbative result 
than its Standard Perturbation Theory (SPT) counterpart. Unfortunately, there are 
several different prescriptions which application is, in some cases, ambiguous. In this 
work, we will calculate the renormalization constants for the B-parameter with the 
following definitions of the coupling constant: 



^In eq.(19) of ref.p^ there is a misprint: the correct value of 722^ is 722^ — —8/3. 
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SPT: Standard Perturbation Theory, af^^ = 5^/(4-71) = 6/(4 tt/?). It is well known 
that this coupling leads to poor convergent perturbative series due to tadpole 
effects We will use this coupling only for the sake of comparison. 

NBPT(n): Naive Boosted Perturbation Theory based on the plaquette; as = 
af^-^/tiQ, where uq = {l/3Tr{Un))^^'^ , was first introduced by Parisi In our 
simulation at /? = 6.0, l/3Tr{Un) = 0.5937. 

NBPT(kc): Naive Boosted Perturbation Theory based on the the critical Wilson 
hopping parameter, Kc', a^^^'^^^"^ = u^^^/uq, where uq = {8kc)~^- In our 
simulation at /? = 6.0, Kc = 0.14543 [||]. 

BPT(ln(n)): Non~perturbative Boosted Perturbation Theory based on the plaque- 
tte's logarithm, ln(n). Following ref.|19|, to obtain the coupling af^'^^^^^°^\q*) 



1. we solve for ay(3.41/a) the expansion of the plaquette's logarithm 

1 Att 
ln{-Tr{Un)) = — -ay(3.41/a) [1 - (1.19 + 0.025 nj) 0^(3.41/0) ] ; 
o o 

2. then we evolve ay(3.41/a) down to the Lepage-Mackenzie's optimal scale 
q* (see below) through the two-loop renormalization group equation with 
the number of flavours nj = because our simulation is performed in the 
quenched approximation. 

Having defined the coupling, it remains to fix the scale at which it will be evalu- 
ated. Since af^"^, a^^^"^^°'* and a^^^'^^'^"^ do not run, no scale setting is necessary 
in contrast to the Lepage-Mackenzie's af^'^^^'^^°^\q*). Taking into account that q* is 
only meant to be a typical scale of the process under consideration, we may guess the 
scale through some physical arguments or, when the two-loop contributions are known, 
choose q* so that the one-loop coefficient of the perturbative expansion vanishes. How- 
ever, Lepage and Mackenzie have suggested a prescription to set q* which consists in 
calculating the expectation value of In(g^) in the one-loop perturbative lattice contribu- 
tion (the integrand of the integral defining the coefficient in front of g'^ /IQtt'^ in eqs.(pO[) 



to (|23D and eqs.(|26D to (|29| ) ). It is clear that for divergent dependent) renormal- 
ization constants, q* depends on /i. Hernandez and Hill |2^, however, have estimated 
q* for heavy-light two-fermion operators using the Lepage-Mackenzie's prescription 
without including the divergent terms proportional to ln(;U^a^). Thus q* is indepen- 
dent of the renormalization scale fi. We will use this scheme to determine the optimal 
scale q* for Za- Further, the Lepage-Mackenzie's prescription is non ambiguous only 
if the operator under consideration does not mix with others under renormalization. 
In fact, if there is mixing, one can either calculate all renormalization constants at the 
same scale q* or determine a separate scale for each mixing coefficient |20]. Since, in 



our case, the operators are divergent and do mix under renormalization, we do not 
know how to implement the Lepage-Mackenzie's prescription. In view of these diffi- 
culties, we choose to determine q* for the heavy-light axial current renormalization 
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constant, Za-, following ref.[^, but for four-fermion operators we will vary (f in the 
range 1 < a(f < n, studying the dependence of Bb on this scale. In this study, for 
simplicity, we will use the same scale q* for all mixing coefficients. Moreover, we will 
use the results from the range of q* to estimate a systematic error. 

Following Lepage and Mackenzie, we now construct tadpole improved lattice oper- 
ators by non-perturbatively redefining the hopping parameter k, k = kuq, where the 
mean-field parameter uq can be taken to be the plaquette or 1/8kc. This redefinition 
leads to the rescaling of the standard normalization of lattice quark fields, \/2k, with 
^/u(), \/2K-y/Mo, which is expected to remove tadpole contributions in the fermionic 
sector Ijl^. Notice that tadpole contributions appear even when a good expansion 
parameter is used. To be consistent, the one-loop perturbative tadpole contributions 
must be subtracted from our perturbative expressions by multiplying by I/^/uq ex- 
panded up to one-loop. Moreover, all links operators U{x) that appear in our lattice 
composite operators, should be rescaled with uq, U{x)/uo, for the same reason as 
above. This method, usually referred to as Tadpole Improvement (TD), is expected to 
reorganize the perturbation series in such a way that TD lattice operators have smaller 
discretization errors and the lattice-to-continuum renormalization factors are nearer 
to unity. For the Wilson action, TD has been studied in ref.||l^. Its HQET counterpart 



is described in refs.|24, ^ to which we refer the reader for details. An important point 
to note is that for HQET quarks TD leads to a final result which is equivalent to the 
reduction of the heavy-quark self-energy constant e, already performed in the previ- 
ous section. Therefore, heavy quark fields are already tadpole improved in our case. 
Light quark fields, however, are not. As explained before, TD involves multiplying all 

light-quark fields by y^/y^?/^ where Uq^^ is the one-loop perturbative estimate of 
the mean-field parameter uq. For example, consider a generic heavy-light four-quark 
operator O. Its expression in terms of lattice operators is 



O = (no) 



where, to one-loop order, Uq^'* = 1 -|- D^q- For two-quark heavy-light operators, 

like the Axial current A^, (uq) in eq.(p9|) should be replaced with ^/uo and Dug by 
Dug/2, because now O contains only one light field. We will implement TD for both 
two and four-quark heavy-light operators and for all choices of the coupling constant 
except for SPT. 



As can be seen from eq.(39), the value of q* for the tadpole improved Za depends on 



the choice of the mean-field parameter uq. We have used two definitions of uq, namely 



the plaquette and 1/8kc- Following ref.|25], where the calculation for the Wilson action 
was performed, and taking into account the contribution from the Clover action, we 
readily obtain for r = 1.0, qA^{a) a = 2.29 for = l/3Tr{Un) and q\ (kc) a = 2.63 
for no = 1/8kc, to be compared to the Wilson result g^^(Kc) a = 2.18 [p5| . 

Finally, we would like to discuss our estimation of the systematic error due to the 
truncation of the perturbative series. The matching factors Z's are the products of the 
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renormalization factors calculated in steps 1 to 3 in the Introduction, which have been 
computed independently up to 0{as). In fact, the QCD operator Oll at the scale rrib 
can be written in terms of HQET operators as follows 
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where ^ is the renormalization scale in the HQET and the functions Xi[^) and Xs{^) 
have been calculated at NLO in refs.[p6| (see also p^, The HQET operators Oll 



and Oj^j^ can, in turn, be expressed in terms of lattice HQET operators using eqs.(20) 
and (|^), 
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Adding all together, we get the QCD-HQET lattice matching for Oll-, 
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(42) 



with obvious notation. We can organize the renormalization constants in two ways: 
including O(a^) contributions coming from the products of the factors expanded sepa- 
rately to 0{as) (method Mi) or excluding O(a^) terms in eq.(|4^) (method M2). Notice 
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that the method Mi, in contrast to method M2, leads to results which depend on the 
renormalization scheme used in the definition of the intermediate continuum HQET 
operators. Since the lattice perturbative corrections are large, different prescriptions 
may result in rather big discrepancies in the final value of Bs^rnb)- We will view these 
differences as an estimation of the systematic error on Bb due to unknown higher 
order contributions to the perturbative renormalization constants. In refs.p^, p^ ], the 
method Mi excluded terms of O(a^) coming from the operator Ofj^, i.e. they take 
Yf^ = Yll = Yrr = Yp = in eg. (|4^) . This is inconsistent because there is no reason 
for excluding these contributions when estimating the systematic error. Numerically 
the difference is very small except for the contribution proportional to Op* due to the 
large value of its matrix element between B-meson states. This has been measured on 
the lattice using the same sample of gauge configurations as for Oll-, Orr, On and 



0£^. We refer the reader to ref.|27] for details. 

Now we give our final results. In Tables 6 and 7, we present the numerical values 
of the renormalization constants and the B-parameter BB{mh) for different choices of 
the coupling constant, the Lepage-Mackenzie's scale q* and different ways in which we 
can organize the perturbative series (methods Mi and M2) Our results have been 
obtained using as(//) at NLO, = 5 GeV, ~ ^00 MeV and four active quark 

flavours, nj = 4. As can be seen, Zop reduces the value of Bb from method M2 by 
an amount smaller than the statistical errors (3-4%) and resulting in a reduction of 
the discrepancy between methods Mi and M2. Further, the inclusion of the correct 
Dji gives a value of which is roughly one-half of the one used in refs.[p^, |2^ . 
Its non-negligible effect on Bb is to increase its value by an amount slightly larger 
than the statistical errors (about 6-8 %). By comparing the results obtained without 
and with TD in options 2 to 5 of Tables 6 and 7 for NBPT couplings, it is clear that 
the Tadpole improvement of the operators significantly reduces the difference between 
methods Mi and M2 and hence the systematic error. In fact, when TD is implemented 
results from methods Mi and M2 are completely compatible within statistical errors. 
We have also implemented TD for a^^^^^^^°^\q*) with two definitions of the mean- 
field parameter uq: using the plaquette or the hopping parameter. Again we found that 
the two definitions give results which are in perfect agreement within statistical errors. 
Moreover, the values of Bb obtained with a^^'^^^'^^'^^^q*) are almost independent of q* 
in the range 2 < q* < n and compatible within statistical errors with those calculated 
with a^^P'^i^) and a^^^'^^'^'^) for both methods Mi and M2. For this reason, we 
estimate our final value of the S-parameter for the coupling a^^'^^^'^^^)) from the 
results in the interval 2 < g* < vr. Finally, we have made the exercise of removing TD 
for this coupling. The main consequences are the following: 

1. the value of the B parameter significantly depends on q* . For instance, from 
method Mi, Bb varies from 0.74(5) at q*a = 2 to 0.86(5) at q*a = vr. 

2. the discrepancy between methods Mi and M2 is much larger than with TD. 
For instance, from method M2, Bb = 0.65(5) at q*a = 2 and Bb = 0.80(5) at 
q*a = vr, to be compared to the values from method Mi given above. 
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Options 


Zol 


Zor 


ZOn 


Zos 


Zop 


Za 






Ml 


0.7935 


-0.0196 


-0.0812 


-0.1130 


-0.0013 


0.9045 


0.0796 


Ms 


0.7658 


-0.0236 


-0.1008 


-0.1229 


0.0000 


0.8999 


NBPT(D) 

as ^ ' 




0.7022 


-0.0330 


-0.13(37 


-0.10G3 


-0.0023 


0.8145 


0.1341 


M2 


0.6539 


-0.0397 


-0.1698 


-0.1229 


0.0000 


0.8067 


NBPTUc) 

as 


Ml 


0.6824 


-0.0359 


-0.1487 


-0.1048 


-0.0025 


0.7950 


0.1458 


Ms 


0.6297 


-0.0432 


-0.1847 


-0.1229 


0.0000 


0.7866 


^^BPnu) TD 


Ml 


0.7302 


-0.0330 


-0.1367 


-0.1084 


-0.0023 


0.8312 


0.1341 


M2 


0.7163 


-0.0397 


-0.1698 


-0.1229 


0.0000 


0.8264 


^NBPn.^) TD 


Ml 


0.7446 


-0.0359 


-0.1487 


-0.1111 


-0.0025 


0.8326 


0.1458 


M2 


0.7387 


-0.0432 


-0.1847 


-0.1229 


0.0000 


0.8282 


* 

q a 


BPT(\n(p)) 

as ^ " 




1.000 


0.2449 


Ml 


0.7092 


-0.0603 


-0.2497 


-0.1148 


-0.0041 


0.7993 


Ms 


0.6977 


-0.0725 


-0.3102 


-0.1229 


0.0000 


0.7957 


2.000 


0.1812 


Ml 


0.7061 


-0.0396 


-0.1636 


-0.1111 


-0.0030 


0.7993 


Ms 


0.7099 


-0.0480 


-0.2054 


-0.1184 


0.0000 


0.7957 


2.290 


0.1727 


Ml 


0.7026 


-0.0370 


-0.1528 


-0.1104 


-0.0028 


0.7993 


Ms 


0.7083 


-0.0449 


-0.1921 


-0.1177 


0.0000 


0.7957 


2.630 


0.1648 


Ml 


0.6986 


-0.0347 


-0.1430 


-0.1098 


-0.0027 


0.7993 


Ms 


0.7060 


-0.0421 


-0.1800 


-0.1170 


0.0000 


0.7957 


3.000 


0.1580 


Ml 


0.6944 


-0.0327 


-0.1347 


-0.1093 


-0.0025 


0.7993 


Ms 


0.7033 


-0.0397 


-0.1698 


-0.1163 


0.0000 


0.7957 


TT 


0.1557 


Ml 


0.0929 


-0.0320 


-0.1320 


-0.1091 


-0.0025 


0.7993 


Ms 


0.7023 


-0.0389 


-0.1664 


-0.1161 


0.0000 


0.7957 



Table 6: Renormalization constants for different choices of the lattice coupling constants, q* 
and options, see the text. The value of Za for the coupling af^^^^^^°^^ has been calculated at 
— 2.63. For this coupling, TD is computed using Uq — I/Skc- 

3. the average value of Bb is smaller than the TD one. 

Prom Table 7, our best estimates of BB{mb) and the renormalization invariant 
S-parameter, Bb, are 

BB{mb) = 0.81 ± 0.05 ± 0.03 ± 0.02 

BB{mb) = 1.29 ± 0.08 ± 0.05 ± 0.03 (43) 

where the first error is statistical, the second is an estimate of the uncertainty due to 
the choice of the coupling constant and the optimal scale q* and the third takes into 
account the contribution of higher-order terms to the renormalization constants. We 
have computed i?B(m(,) using the NLO formulae. 
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Options 










Bop 




0.0796 


Ml 


0.912(43) 


-0.023(1) 


-0.097(6) 


0.083(4) 


-0.013(0) 


0.862(44) 


M2 


0.889(42) 


-0.027(1) 


-0.122(7) 


0.091(4) 


0.000(0) 


0.831(43) 


NBPTin) 
as 

0.1341 


Ml 


0.995(47) 


-0.047(2) 


-0.202(13) 


0.096(5) 


-0.027(0) 


0.816(49) 


M2 


0.945(45) 


-0.057(2) 


-0.256(16) 


0.113(6) 


0.000(0) 


0.745(48) 


NBPT(kc) 

as 
0.1458 


Ml 


1.015(48) 


-0.053(2) 


-0.231(15) 


0.099(5) 


-0.030(0) 


0.800(51) 


M2 


0.957(45) 


-0.066(3) 


-0.293(19) 


0.119(6) 


0.000(0) 


0.718(50) 


NBPT(D) 

as TD 
0.1341 


Ml 


0.993(47) 


-0.045(2) 


-0.194(12) 


0.094(5) 


-0.026(0) 


0.823(49) 


M2 


0.986(47) 


-0.055(2) 


-0.244(15) 


0.108(5) 


0.000(0) 


0.796(50) 


as ID 
0.1458 


Ml 


1.010(48) 


-0.049(2) 


-0.210(13) 


0.096(5) 


-0.028(0) 


0.819(50) 


M2 


1.012(48) 


-0.059(2) 


-0.264(17) 


0.107(5) 


0.000(0) 


0.797(51) 


* 

q a 


^BPT(^nin)^ 




1.0000 


0.2449 


Ml 


1.043(49) 


-0.089( 4) 


-0.383(25) 


0.108( 5) 


-0.050( 0) 


0.629(56) 


M2 


1.036(49) 


-0.108( 5) 


-0.480(31) 


0.116( 6) 


0.000( 0) 


0.564(59) 


2.0000 


0.1812 


Ml 


1.039(49) 


-0.058( 2) 


-0.251(16) 


0.104( 5) 


-0.036( 0) 


0.798(52) 


M2 


1.054(50) 


-0.071( 3) 


-0.318(20) 


0.112( 5) 


0.000( 0) 


0.777(55) 


2.2900 


0.1727 


Ml 


1.034(49) 


-0.054( 2) 


-0.234(15) 


0.104( 5) 


-0.034( 0) 


0.814(52) 


M2 


1.052(50) 


-0.067( 3) 


-0.297(19) 


0.112( 5) 


0.000( 0) 


0.799(54) 


2.6300 


0.1648 


Ml 


1.028(49) 


-0.051( 2) 


-0.219(14) 


0.103( 5) 


-0.032( 0) 


0.828(51) 


M2 


1.048(50) 


-0.062( 2) 


-0.279(18) 


0.111( 5) 


0.000( 0) 


0.818(53) 


3.0000 


0.1580 


Ml 


1.022(48) 


-0.048( 2) 


-0.207(13) 


0.103( 5) 


-0.031( 0) 


0.839(51) 


M2 


1.044(50) 


-0.059( 2) 


-0.263(17) 


0.110( 5) 


0.000( 0) 


0.833(53) 


vr 


0.1557 


Ml 


1.019(48) 


-0.047( 2) 


-0.202(13) 


0.102( 5) 


-0.030( 0) 


0.842(50) 


M2 


1.043(49) 


-0.058( 2) 


-0.258(16) 


0.110( 5) 


0.000( 0) 


0.837(53) 



Table 7: Values of each operator contribution to the B -parameter, and the total result, for 
different choices of the lattice coupling constants, q* and options, see the text. 



We have repeated the same analysis with the lattice data for the B-B mixing from 
the UKQCD Collaboration obtained using the SW-Clover and HQET lattice actions 
at (3 = 6.2 [^]. The conclusions presented above remain the same and the final results 
turn out to be in perfect agreement with ours in eq.(^), 

BBinib) = 0.79 ± 0.04 ± 0.03 ± 0.02 

BBinib) = 1.26 ± 0.06 ± 0.05 ± 0.03 (44) 

Note also that we do not see any dependence of our results on the lattice spacing a, 
within errors. 



These numbers are to be compared with our previous result |22|: Bsinib) = 
1.08 lb 0.06 lb 0.08, obtained with the same lattice data but with the incorrect values of 
the renormalization constants and without tadpole improvement for the four-fermion 
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operators. Unlike our previous result, our best estimates in eq.(|4^) are now in good 
agreement with previous determinations of the 5-parameter calculated by extrapolat- 
ing Wilson data from the charm to the bottom mass (see Table 4 of ref.[22|). 



8 Conclusions 

In this paper we have determine the expressions of all AB = 2 and AB = four- 
fermion HQET continuum operators in terms of HQET lattice ones using perturbation 
theory up to one loop. This calculation is a necessary ingredient in the procedure to 
measure the unknown values of some important matrix elements of QCD four-fermion 
operators by simulating them with lattice HQET. The calculation has been performed 
in four continuum renormalization schemes, NDR, DRED I, DRED II and HV, and for 
three lattice actions, Wilson, SW-Clover with rotated fermion fields and SW-Clover 
with rotated propagators. 

We have also discussed the subtleties associated with the renormalization scheme 
dependence of our results. The effects of evanescent operators which appear in inter- 
mediate steps of the computation have been identified and clarified. 

As a byproduct of our calculation, we have also obtained the one-loop anomalous 
dimension matrices for both the AB = 2 and AB = HQET operators. 

We have found and corrected an error in the one-loop result quoted in |jl2| for the 
lattice-continuum matching of the operator Oll, which matrix element between B- 
meson states determines the value of the Bb parameter of the B-B mixing. Then we 



have reanalyzed the lattice data measured in ref.[22| with the correct renormalization 
constant. We have used Boosted perturbation theory with tadpole improvement to 
reduce the systematic uncertainties in the renormalization constants due to the trun- 
cation of the perturbative series. Our value for the Bb parameter turns out to be 
larger and has a smaller systematic error than in previous studies with static heavy 
quarks. Further, it is compatible with the results obtained by extrapolating lattice 
Wilson data. 
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Appendix A Analytical expressions of the con- 
stants. 

In this appendix we give the analytical expressions for the lattice constants which enter 
the continuum-lattice matching relations. 

= 2 operators. 

Operator Oll 



Dll = -l/3v - 10/3di - l/3c - 4/3e - 4/3/ 

= - 1/3 + 8/3 (Z + m) + 1/3 s - 4/3/^ - 10/3 n 

Dii = h 

Drr = 4/3 ^« 

£>^^ = 4/3u;^ - 2/3(/ + m) + l/3s 

Drr = Sr 

Dn = 2d2 

= -4d^ - 2q + 2h (A.l) 

Operator Olr 



Dlr = + 1/6 Ji - 7/3(ii + l/6c- 4/3e - 4/3/ 

Diji = + 13/6 +m) + 1/4 s + 1/6 J2 -4/3/^ - 7/3 n 

Di'j, = 3/4 5^ + 1/6 J3 

Df^L = 2/3w + Ji + 2di + c 

Dfii = 2/3u;^ - l/3(/ + m) + l/6s + J2 + 2n 

5|f = l/25R + h 

Dm = d2 

Dij = -2/ - q + h (A.2) 

Operator 



= -2/9u + 8/9 Ji - 4/3(ii + 2/3c - 4/3e - 4/3/ 

Dll = -2/9u^ + 8/3(Z + m) + 2/9s + 8/9 J2 - 4/3/^ - 4/3n 

Dll' = 2/3 -^L + 8/9 J3 

Dll = l/36u + 2/9 Ji + l/2di + l/4c 

DIl = 1/36 - 1/36 s + 2/9 J2 + 1/2 n 
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-1/12 5l + 2/9 J3 
-l/3tt; 

-l/'iw^ + l/Q{l + m) - l/12s 
-l/4<5i? 
-1/4^2 

l/2d^ + l/4g - l/4/i 
Operator O;^^^ 

w + 1/6 Ji - 7/3 di + 1/6 c- 4/3 e - 4/3/ 

+ 13/6 (/ + m) + l/4s + I/6J2-4/3/ - 
3/4 Jfi + 1/6 J3 

1/6 w; + 1/4 Ji + 1/2 di + 1/4 c 
1/6k;^ - 1/12 (Z + m) + l/24s + 1/4 J2 + 1/2 
1/8 + 1/4 J3 

-l/4d2 

l/2d^ + l/Aq - 1/Ah 

AB = operators. 

Operator Qll 

Ell = -8/3di - 4/3e - 4/3/ 

eIl = 8/3(/ + m) - 8/3n - 4/3/^ 

EtLL = Ji + 2di + c 

EtiL = J2 + 2n 

Eti^L = J3 

e4l = -^w 

Et^i = -Aw^ + 2{l + m) - s 

Et^RL = - 3 
Er = -4/3d2 
Eji = 4/3 (2d^ + q - h) 

EtR = -d2 

Et^ji = 2d^ + q - h 

Operator Qlr 



Drr 

D'rr 

^rr 
Dp 



^lr 

r,si 
^lr 



D 



SII 
LR 



Drl 



D 



D 



RL 
II 
RL 



Dq 
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Elr = -8/3di - 4/3e - 4/3/ 

Ei^ = 8/3(/ + m) - 8/3n - 4/3/ 

EtLR = V + 2di + c 

Et[ji = - s + 2n 

Eti'^ = -3Sl 

EtRL = 

EtRL = -'^w^ + 2{l + m) - s 

EtRL = - 3 <5r 

Es = -i/3d2 

E^S = 4/3 (2d^ + q - h) 

Ets = -d2 

Et^S = 2d^ + q - h (A.6) 



Operator Qfj^ 



Ell 




-8/3di - 4/3 e - 4/3/ 


^LL 




8/3(/ + m) - 4/3/^ - 8/3 n 


EtlL 




4/3 Ji - l/3v + 2di + c 


Etfl 




-l/Svi + l/3s + 4/3 J2 + 2n 


EtlU 




Sl + 4/3 J3 


EtLR 




-1/3 Ji + 1/3 1; 


Et^LR 




1/3 - 1/3 s - 1/3 J2 


EtLR 




-Sl - 1/3 J3 


EtRL 




— w 


EtRL 




-w^ + 1/2 (Z + m) - l/4s 


EtRL 




-3/^Sr 


Et 




-4/3d2 


Ejt 




4/3 (2(i^ + q - h) 


EtT 




-d2 


Etip 




2d^ + q - h 



Operator Q^r 



= -8/3di - 4/3e - 4/3/ 
28 



Constant 


r = 1 


r = 0.75 


r = 0.50 


r = 0.25 


r = 


Jl 


-4.85 


-4.93 


-5.18 


-6.08 


-8.24 


J2 


-0.16 


-0.21 


-0.28 


-0.25 


0.00 


J3 


-0.11 


-0.09 


-0.06 


-0.01 


0.00 


h 


0.20 


0.15 


0.09 


0.02 


0.00 


6r 


-1.23 


-0.58 


-0.17 


-0.01 


0.00 



Table A.l: Numerical values of the new integrals for various values of the Wilson parameter 
r. 



s 



E^LR 

EtRR 



Et 



Et 



RR 
II 
RR 

Eu 

E: 



u 



Etu 
Et\j 



8/3 (Z + m) - 8/3n - 4/3/^ 
Jl + 2di + c 
J2 + 2n 

h 

— w 

-w^ + l/2(/ + m) - l/4s 

-3/4 5h 
-4/3d2 

4/3(2/ + q - h) 
2d^ + q - h 



(A. 



where the expressions and values for the constants di, d2, f, v, w, c and e which 
determine the matching in the Wilson case, can be found in ref. |21|. The Clover con- 
stants I, m, n, q, h, s, d^ , f^, and are defined in ref.[12| whereas the analytical 



expressions of 6l and 6r are given in ref.[13|. We have found small numerical discrep- 
ancies between our estimates for d^ , f^, h and and the values quoted in ref. . 
Our results for r = 1.0 are d^ = -4.13, = -3.64, h = -9.97 and = -6.72 to 
be compared to d^ = —4.04, = —3.63, h = —9.88 and = —6.69 from Table 3 
of ref.p^]. We have used the former to obtain all our results which agree with those 



recently and independently obtained in ref.|29]. 

The only new constants are Ji , J2 and J3 defined as 



Jl 
J2 





/ d^k 

J—7T 


"0(1 -A;2) 




k" 




/ d^k 


" r^Al 




8A1A2 



A4 



+ 



4AiAi 4AiAi 
^^(A4- A5) ' 



r^Af 
A2 



2A2 
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J?, 



where 



TT 



16Ai 



Ai 
A2 
A4 
As 



4A2 



(Ai 



2 9^ 
2 



^ sin^g^ + A? 



^ sin^ 
^ sin^ sm 



2 9m 
2 



(A.IO) 



In Table A.l, we give the numerical values of Jj, i = 1,2,3 for several values of the 
Wilson parameter r. We also give the constants 5l and because in ref.|13| they are 
not presented explicitly. 
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Figure 1: Feynman Diagrams for the continuum-lattice matching. Up- down replicas and 
t -channel diagrams are not shown. For the AB = operators only s -channel diagrams 
contribute, penguin diagrams are not considered and the arrows in the lower fermionic lines 
must be inverted. 
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